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CHAPTER I 


INTRODUCTION 

1.1 Motivation and Re levan ce o f Essay 

The balloon system has been employed extensively in the 
past as a means of conducting research in the earth's 
atmosphere. The main disadvantage in the use of balloon 
systems occurs in those experiments where it is necessary to 
either stabilize the observation platform or use some observer 
system to perdict its attitude as a function of time. 
Stabilization or prediction is necessary at times in order 
to process the experimental data collected by various research 
instruments which are mo\mted on the observational platform. 

Platform stabilization can be accomplished by means of 
control systems. However, these systems are usually complex 
and result in increase cost and additional platform weight 
which reduces platform payload. Another method is to allow 
the platform to oscillate freely in space and predict its 
attitude. Thus, it is important that designers of balloon 
systems have a understanding of the motion of balloon borne 
observation platforms in order to determine what kinds of 
auxilliary systems are necessary for accurate and economical 
data collection- 

The motion of balloon borne observation platforms during 
ascent and decent have been discussed in several papers 
(Ref. 1,2). In this essay the effect of wind ousts which 
result in forces acting externally on the ballooon system will 
be studied. These forces affect the nature of the motion of 


the observation platfoi'm while the balloon is at float 
altitude. 

^ ^ Objective o f Essay 

The objective of this essay is to deter?iiine the effect 
of wind gusts on the magnitude of the pendulation angles of 
a balloon borne observation platform. A system mathematical 
model will be developed and the solution of this model (in 
conduction with data gathered by NASA during a flight mission) 
will be used to determine the magnitude of the observation 
platform pendulation angles. 
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CHAPTER II 

DEVELOPMENT OF MATHEMATICAL MODEL 
2 . 1 LACATE Experiment 

The National Aeronautics and Space Administration 
conducted a high altitude balloon experiment called LACATE 
(Lower Atmosphere Composition and Temperature Experiment) . 
This experiment employed an infrared radiometer to sense 
remotely , vertical profiles of concentrations of selected 
atmospheric trace constituents and temperature. 

The balloon system used for the mission is shown in 
Pig. (1) . The system consists of the following: 

a. A 39 million cubic feet zero pressure balloon. 

b. Radar reflector. 

c. Recovery parachute. 

d. Flight control electronics. 

e. Radiometer. 

f. A platform containing the research payload. 

When the balloon attained float altitude, data v/as 

gathered by the various instruments and telemetered to 
ground control. The radiometer line of sight was scanned 
vertically across the horizon at approximately 0.25^ 
per second, requiring 30 seconds to acquire a complete 
radiance profile. At the end of the mission the platform 
system was separated from the balloon and returned to earth 
by means of the parachute. 

2 - 2 Idealization of Balloon Platform System 


The actual motion of the balloon system once it reaches 
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float altitude is very complex and involves various types 
of oscillations, including bounce {vertical oscillations), 
pendulations (in plane motion) , spin (rotation) , and horizontal 
translation . 

In previous work (Ref. 3) the LACATE system was idealized 
as shown in Fig. (2). Each balloon subsystem was treated as an 
equivalent rigid body. The mass of the entire system was 
lumped at the center of gravity of the balloon and positions 
1, 2, and 3 as shown in Fig. (2), Euler angles can be choosen to 
measure spin and pendulation in two mutually perpendicular 
planes. It was shown in (Ref. 4) that by choosing the proper 
set of Euler angles and assuming small displacements the 
pendulation motion uncouples in two mutually perpendicular 
planes, thus simplifying the form of the mathematical model. 
Details of that development can be found in (Ref. *4) , 

Since the pendulation motion uncouples, the balloon 
platform system for this study will be further idealized 
as shown in Fig. (3) in order to develope the form of the 
system mathematical model. This model in the mutually 
perpenducular plane (i.e., y-z plane) is the same, with 

®1 “ ~ 

For purpose of this study the following assumptions 
are made : 

1. The distributed balloon subsystem will be lumped 
into two subsystems and treated as equivilant 
rigid particles. 

2. The cables will be treated as inflexible and 


inextensible 
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3. The altitude of the support point 0 will assumed 
to be constant during the entire period of 
observation; i.e. , y = 0. 

4 . The moments of interia of each subsystem about 
their center of gravity will be neglected. 

5. Viscous drag forces, viscous drag torque and 
support friction will be neglected. 

This idealization enables one to treat the balloon system 
as a double pendulum with a moving (accelerating) support. 

The acceleration of the support is due to wind gust forces 
which act on the balloon. These accelerations cause equiva- 
lent excitation forces that affect system response. The 
main purpose of this work is to study the effect of these 
forces on the angles 0-]^, 

2 . 3 Gene ra lized Coordinates 

The generalized coordinates for a given system are those 
coordinates which are employed to specify the configuration of 
the system at any instant of time. In any mechanical system 
the number of degrees of freedom of the system coincides with 
the minimum number of independent coordinates necessay to 
describe the system uniquely. In the case of the idealized 
planar lumped parameter system shown in Fig. (3), two 
generalized coordinates, angles (0^ and 02 or and 
necessary to specify the configuration. For this study 
generalized coordinates will be employed in order to facil- 
take the use of Lagrange's Equation for developing the 
mathematical model. 


2.4 
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■uagrangc 


JKquation 


^he diftersntial equations governing the motion of the 

syste. Shown in Pi,. ,3, win be developed by employing 

Lagrange's equation. The form of fho 

■ine rorm of the equation used in this 

study is given as follows; i.e.. 


where 


_d ^ 

d t 5 i 

dgk 


9qb 


= Q 


k' 


k=l,2, . . 


. n 


( 2 - 1 ) 


L = T - V = Lagrangian, ^2_2) 

T = kinetic energy of the system, 

V = potential energy of the system, 

^k ~ 9®riS3rali2ed coordinate, 

^k ~ Sfsneralized velocity, 

number of generalized coordinates, and 
Ok = the nonconservative generalized forces. 

represents a set of n simultaneous differential 
equations which describe the motion of a holonomio system. 

In this study support friction and drag forces are 
neglected therefore = q, 

/V 

Ihe two main advantages for using Lagrange’s equation are- 
1- The internal reaction forces do no work during the 
■notion and therefore can be neglected. 

The energy terms can be computed in a straight- 
forward manner. 
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2 . 5 System Lagrangian 


The kinetic energy of the system shown in B'ig. (3) 
(assuming small angles) is given as follows; i.e., 

T = ism^(x + + 1^0^ + 1262 )^' (2-3) 

where 


T = kinetic energy of complete system, 
m^ = mass of subsystem 1 , 
m 2 = mass of subsystem 2 , 

1^ = distance from support point O to mass 

I 2 = distance from mass m^|^ to mass m 2 , 

X = translation velocity of support point O, 

6 ^^ = angular velocity, and 

62 “ angular velocity. 

The potential energy of the system is given as follows 

V - m^^glj^ (l-cosG^) + m 2 g (Ij^ (1-cos 0]^) 

+ I 2 (i^-cos02) ) (2-4) 

where 

V = potential energy fo the complete system, and 
g = acceleration of gravity. 

Substitution of Eqs. (2-3) and (2-4) in-Tri Eq. (2-2), 
expanding and simplifing terms yields the system Ijagrangian 
and is given as follows; i.e.. 
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L = 


’ 2 


+ 2L^XGj^ + + Hm2 (X^ + 2 XLjl0^ 

+2XL202 + ^^i^2®1®2 + ^2^2^ ■■ 

+ m^gLj^cosO^ - m^gL^ + in2gLj^cos0^ - ni^gL^ 


+ m gL_cos( 
2 2 


{2-5) 


2 • 6 Equations of Motion 

The differential equations governing the motion of the 
system shown in Fig. (3) are obtained by substituting 
Eq. (2-5) into Eq. (2-1) with q = 0 and q = ^ . The 

results after rearranging terms and assuming small angles 
are given as follows; i.e., 


(m^L^ + 


2i 

•2-1- 


m„Lf ) 0. 






~ ^^1^1 ^ ^2^1^^' 

( 2 - 6 ) 

^"' 2 ^ 1 ^ 2^®1 ■*■ "' 2 ^^ 2°2 " “ "’ 2 ^ 2 ^- 

Eqs, (2-6) and (2-7) are the differential equations 
governing the motion of the idealized balloon system. The 


equations of motion in the mutually perpendicular plane 


are the same, however 0^ = and 0^ = O.Kse equations 
can be written as: 


>2 2 " 

(miLi + m2L^)i^^ + m2Lj^L2i|;2 + + m2gL^)i|j^ = 

- + ro 2 L^)X, and 


( 2 - 8 ) 
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(m2Lj,L2H'j^ + (m2L2)i^2 ” “*^2^2^' (2-9) 

One method of solving these equations will be discussed 
in the following chaptoi*. 


13 


CHAPTER III 

RESPONSE OF BALLOON BORN OBSERVATION PLATFORM 

3.1 Development of Lumped Parameter Modal Equations 

The method of modal analysis can be used to transform 
the simultaneous coupled differential equations of motion 
of a lumped system into a set of uncoupled differential 
equations, (Ref. 5). These resulting equations can easily 
be solved to obtain the response as a function of various 
initial conditions and excitations. 

The mathematical model for any generally linear lumped 
parcimeter mechanical system without damping can be written 
in matrix form as follows; i.e., 

Mq + Kq = F, (3-1) 

where 


M = mass matrix, 

K = stiffnes matrix, 

q = vector of generalized coordinates, and 
F = forcing function. 

To use the method of modal analysis, it is necessay 
to solve the eigenvalue 'problem associated with thfe ' 
homogeneous system described by Eq. (3-1) . The eigenvalue 
problem can be expressed as follows; i.e. , 

w?Mu^ = Ku^ 


(3-2". 
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where 


2 th 

0)^ = i eigenvalue , and 


. th 

= 1 eigenvector . 


The eigenvectors 'can be normalized such that: 


uTmu . = 6 . . 

1 3 13 ^ 

S^Ku, = 


( 3 - 3 ) 

( 3 - 4 ) 


where 


th 


- i normalized eigenvector, 


~ . th 

Uj = 3 normalized eigenvector, and 


= Kronecker Delta. 


= 1 i=j 

= 0 iyij 


It can be shown that: 


u. = C. li. , 


(. 3 - 5 ) 


where 


2 T 

Cf = u.Mu. . 
1 11 


( 3 - 6 ) 


The resulting modal matrix U is such that; 


U MU = I , 


and 


( 3 - 7 ) 


T 2 

U^KU = u) , 
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(3-8) 

where 

U = [u3^U2l , (3_9) 

1 = identity matrix, and 

2 

w = diagonal matrix of the eigenvalues. 

The non-homogeneous solution of Eq. (3-1) can now be 
described as follows; i.e. , 

^ (3-10) 

where 


n == column matrix consisting of a set of time 
dependent geneialized coordinates, and 
q = column matrix of generalized displacements. 

Substitution of Eq. (3-10) into Eq. (3-1) yields; 

MUii + KUn = F. ^ 2-1 

» 

Premultiplying both sides of Eq. (3-11) by yields 
the following; i.e., 

U^MUii + u\un = U^F (3-1 

Introduction of Eqs. (3-7) and (3-8) into Eq. (3-12) 
gives the following expression; i.e., 

.. , 2 

n + 10 n = N, 


where 


(3-13) 
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T 

N = U F. 


(3-14) 


Eq. (3-13) represents a set of n uncoupled differential 
equations of the form: 


ii^(t) + wfnjL(t) = 


N, 


i = 1 , 2 , 


,n. 


(3-15) 


These equations have the form of the differential equations 
describing the motion of n, undamped, uncoupled single 
degree of freedom systems. 

Eq. (3-15) can be solved by means of the Laplace 
transform method, this yields: 


1i(t) = (t) sinwj. (t-T)dx + ( 0) cosu)j_t 

i = l,2,...,n. ( 3-16) 


0 


+ n, 10) si-pi- 

1 U.'.; 


Assuming zero initial conditions, Eq. (3-16) 
written as follows; i.e.. 


Hi(t) = -JN^(T)sino)^(t-T)dT. 


can be 


(3-17) 


This equation is known as the convolution integral. The 
response q can be determined by introducing Eq. (3-17) 
into Eq. (3-10) 


3 . 2 Response of Idealized Balloon System 

The Eqs. (2-6) , (2-7) , (2-8) , and (2-9) for the system 
shown in Fig. (3) can be written as follows: 
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^ 2 ^ 1^2 


,2 

'’^2^2 


(in^^l 2^"^it\2l ^ ) ^ 

- (m2l2)x 

or in abbreviated form as : 


Mq + Kq = F, 







^1 


mj^gl^+ra2gl2^ 0 


^2 

^2 

4- 

0 ^2^^ 2. 


^2 


(3-18) 


(3-19) 


where 


M 


K 


mil 

"»21 




"*22 


0 

^22 


,and 


- (ra^l^+m2lj^) X 
- (m lo)x 

A non trivial solution of the eigenvalue problem given 
by Eq. (3-2) will result if and only if the determinant of the 
coefficients vanish; i.e., 

2 

|K - to Mj = 0 . 


(3-20) 
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Introduction of the proper elements of Eq. (3-19) 
into Eq. (3-20) yields. 


11 X 


■^12^'^ 


‘ = 0 . 


-m2iw^ 


^22 


(3-21) 


Solving Eq. (3-21) for results in two real roots given 
by the expression: 



(kiimii4-k22inii) 

^™11"'22""'i2”^21^ 


where 


± 1/2 

V 


(kii^22-"^22^1l)^ 

(miim22-mi2in2i)2 


"^^ 11^22 ( 3 - 22 ) 

Ki "' 22 -"’ 12"'21 


0 )^ = natural frequency, and 
i = 1,2. 

The eigenvectors associated with each natural frequency 
can now be written as; i.e., 

1 

""l = 2 ' 

kll -‘^1^11 

a)^mi2 


(3-23) 
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and 


^2 


jjJLVM 

w|mi2 


where 


(3-24) 


- eigenvector associated with and 

^2 “ eigenvector associated with 

The characteristic vectors and U 2 are the modal 
vectors and represent the natural modes of oscillation of 
the system. Normalization of the modal vectors will 
produce the ;normal modes. Computing the values of and 
C 2 according to Eq. (3-6) yields the following,* 



(J'll-wfniii) 

mi /m2 i+mi 2 ) +■ 

ii)imi2 


(kii-wjmii) ^ 


(to?m 


1“'12 


) 


^22 ! (3 



^^11~^2^11^ 

2~ 

“2^^12 


(m2i+mi2) + 


(kii-a)2mii)2 




m 


22 


. (3' 


-25) 


26) 


The resulting normalized modal vectors become: 
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Ui = 


kii - 


“i^l2 


, and 


(3-27) 


“i = 


*^11 ~ 




(3-28) 


Substituting Eqs. (3-27) and (3-28) into Eg. (3-9) 
yields the modal matrix U; i.e., 


U = 


^11 ~‘^l"'ll 




c- 


^11 ~ ^ 2^11 
^2“2"’l2 


(3-29) 


Substitution of the modal matrix U given by Eg. (3-29) 
into Eg. (3-14) yields the follows’ r^g for the components of 
vector N; i.e., 

Nf = - (m|^l^+m 2 l 2 ) xu^j^ - (m 2 l 2 ) XU 2 ^, and’ 

N 2 = - (m]^l^+m2l2) XU2^2 “ (ni2l2^ ^^^22 ' 

where 


Uj^j = components of modal matrix U. 


21 


introducing the the above components into Eq. !3-17) 
yields the following set of expressxons; x.e., 


lift) = C(T)slnu.j^(t't)aT. and 

(3-30) 

0 / t 


= |jjNj(T)sinu.2(t-T)aT. 

(3-31) 


'0 

The system response given by Eg» 


(3-10) can now be 


written follows: 


^ 


and 


(3-32) 


^2 " ^ 1^^^"^21 


(3-33) 


The next 
evaluation of 
of the system 


chapter will discuss the numerical 


and q2 order 


to obtain the response 


shown in Fig. (3). 
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CHAPTER IV 
RESULTS AND SUMMARY 

4 . 1 LACATE Mission Data 

Fig. (1) illustrates the actual LACATE balloon 
system and ^ig. (3) illustrates the system as it was 
idealized for this study. The physical properties of the 
idealized LACATE system are given in Table I. The results 
of the eigenvalue problem from Eq. (3-2) are presented in 
Table II. The resulting natural frequencies, periods, and 
mode shapes are given in Table III. 

4 . 2 Balloon System Horizontal Accelerations 

The flight of the LACATE balloon system was tracked by 
ground radar and the resulting horizontal trajectory is 
shown in Fig. (4). Two typical segments of the balloon's 
flight path were analyzed in this study. Flight segment 
one is taken from the time interval 0£t<^500 seconds, while 
flight segment two is taken from the time interval 
588 9<_t<^638 9 seconds. 

The hoizontal acceleration components and a .2 along 
the balloon's body axis were computed using data from the 
radar tracking station in conjunction with a numerical 
differentiation process, (Ref. 4). Figs. (5) and (6) show 
the body axis acceleration components a^ and ^2 for the 
two flight segments studied. 

4 . 3 Response of Balloon Borne Observation Platform 

The values of and ri 2 fi^om Eqs. (3-30) and (3-31) 


were determined by numerical integration of the convolution 
integral using a trapezoidal rule. (Refer to Appendix I 
for a listing of the computer program.) Fig. (7) -(14) 
give the response of and (with inputs a^^ and a 2 ) for 
each flight segment studied. 

In order to verify the results of the computer program, 
the input a^, from flight segment one was approximated as 
a straight line function for the interval 0 <^t ^200 seconds. 
The resulting approximation is shown in Fig. (15). 

Eq. (3-30) was then integrated analytically with this 
straight line approximation as input for x. The results 
are shown plotted in Fig. (16). Eq. (3-30) was also 
integrated numerically on the computer with the straight 
line function as input for x. (Refer to Appendix II for a 
listing of the computer program.) The results are given 
in Fig. (17) . 

The results from Eqs. (3-32) and (3-33) for 0-j^, 62 ; 
and i |)2 flight segment one are presented in Figs. 

(18) -(23). B_gs. (20) and (23; show the comparision of 
9^, 02 , and tp 2 respectively. The system response 

e^, 02 , and \p 2 for flight segment two are presented in 

Figs. (24)-(29). 

4 • 4 Discussion of Results 

The results in Fig. (4) indicate dramatically the 
presents of wind gusts during the flight of the LACATE 
balloon system. These gusts were sufficient to cause 
significant changes in the direction of the balloon's 


flight path. Moreover, Figs. (5) and (6) indicate that 
because of the wind gusts the magnitude of the acceleration 
a^^ and a 2 varied from -0.123 m/sec^ to 0.095m/sec^ in flight 
segment one and -0.055m/sec to O.lOm/sec in flight segment 
two. 

Figs« (16) and (17) indicate that the results and 
TI 2 obtained from the computer program are in close agree- 
ment with the results obtained from the analytic method. 

This fact enables one to have confidence in the computer 
program and in the numerical integration algorithm. 

Pigs. (7) “(15) show that the response of and TI 2 is 
oscillatory in nature. Moreover, the relative magnitude 
of the oscillation increases or decreases as the input 
acceleration a^^ and a 2 increases or decreases. 

The results in Figs. (20) and (26) indicate that the 
magnitudes of 0^^ and are approximately equal. This is 
also true for and \l >2 shown in Figs. (23) and (29) . 

Figs. (30) -(33) show that the period of oscillation is 
approximately 10 seconds. These facts indicate that the 
idealized balloon system oscillates primarily in the 
most fundamental mode. 

The results in Figs. (18) -(29) indicate that the 
system- does not oscillate about its vertical axis. These 
figures show that the relative offset of oscillation is 
negative when the input is postive and vice versa. The 
figures indicate also, that the relative swing amplitude 
increases and decreases as the input acceleration increases 
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and decreases. 

4 . 5 Summary 

In this study the balloon borne observation platform 
was viewed as a double pendulum with a moving (accelerating) 
support. The acceleration of the support which is due 
to the effect of wind gusts acting on the balloon was 
treated as an idealized source. 

The results of this study indicate that the system 
oscillates primarily in the most fundamental mode. This 
is significant because, the balloon system can therefore 
be treated as a single pendulum with a moving (accelerating) 
support. This will simplify the mathermatical model in 
future idealizations. The results of this study also 
indicate, that due to the effect of wind gusts the 
maximum pendulation angle is of the order 1.0 degree. 
Generally these angles were found to be of the order 0.60 
degrees or less. This is important to designers of future 
balloon systems inorder for them to determine what kinds of 
auxiliary controls and or atitude determination systems 
are necessary for economical and accurate data collection. 



I 
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TABLE I 

Idealized LACATE System Properties 

= 75ft (distance from support point O to mass m^) 

= 15ft (distance from mass m^ to mass m 2 ) 

mi = 1351b (lumped mass) 
m 

m 2 = 3751bj^ (lumped mass) 




TABLE III 


Natural Frequencies, Pe 
Modal Shapes for Ballo 
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3.0624 


2.051 sec 













I I A C AN 


i 

'> 4066 / 




0.0 


100.0 


200.0 300.0 

Time (sec) 

Figure (5) Horizontal Acceleration 
Along Body Axis (Flight 


400.0 


500.0 


of .Support Point O 
Segment One) 




1.5 


1 . 


0 . 


0 




T') 

e 

cr 


fH 

C 


-0.5 

- 1.0 



-1.5 


- 2.0 


-2.5 


-3.0 


0.0 


100.0 200.0 300.0 

Time (sec) 


400.0 


500.0 



Figure (7) Plot of Values (With Input = a^) 


Ul 

K> 


OF i^oo; 











jr* 


<N 

e 

CP 


c 



Values (With Input = aj^) 


O O 

~n 23 
«** 

T3 P> 


o 

^ i-»a. 

O 3D> 

53 i" 


c 

> o 

1“ jj 


4 


k 'a 


Figure (11) Plot of 


tjj 

o\ 


37 


ORJGJWAl P/IGH ?3 
OF POOR QUAUIY 



(12) Plot of rso Values (With Input 



0.75 


0.50 

0.25 

0.00 

-0.25 

-0.50 
-r 

-0.75 
tJ* 

^ - 1.00 

-1.25 

-1.50 
-1.75 
- 2.00 

-2.25 
-2.50 

5889.0 5989.0 6089.0 6189.0 6289.0 6389.0 

Time (sec) 

Figure (13) Plot of n. Values (With Input = a^) 

1 ^ OJ 

CO 



0.004 


0.003 
0.002 
jf, 0.001 

(N 

E 

5 0.000 

OJ 

- 0.001 

- 0.002 

- 0.003 

- 0.004 



(m/ sec 


0.060 


0.040 

0.020 

0.000 

- 0.020 
fN 

-0.040 
^ -0.060 
-0.080 
- 0.100 

- 0.120 
-0.0140 

0.0 100.0 200.0 300.0 400.0 500.0 

Time (sec) 

Figure (15 ) Plot of Straight Line Approximation 
Flight Segment One 



O O 


4 on 


.c. 

o 
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Figure (19) System Response (with Input 








Figure (21) System Response (With Input 
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Figure (22) System Response (With Input = a2) 
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Figure (23) System Response (With Input = a 2 ) 
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Figure (25) System Response (With Input = ) 
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Figure (28 ) System Response (With Input = 32 ) 
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APPENDIX I ^ 

COMPUTER PROGRAM FOR CALCULATING THE FOLLOWING 

Ui, 0)., Uij, 8i, 4 ,^ 


1 . 000 
2 . 000 
• 5 . ooo 

4.000 
0.000 

6 . 000 
7 .. 000 
P.OOO 

9.000 
1 6 . 000 
n . 000 
1 ? . 000 
1 0 . 000 
1 4 .. 000 
1 5 . 000 
16.000 
1 7 . 000 
1 S . 000 
1 q . 000 
20.000 
21 .000 

22.000 
2 '^. 000 
24.000 
25.000 
26.000 
27.000 
2 R. 0(;0 

29.000 

30 . 000 
51 .000 
32.000 
3 '^. 000 

54.000 

55.000 
56.000 

57.000 

ys.ooo 

39.000 C 

40 . 000 

41 .000 

42 . 000 

43.000 

44.000 

45.000 

46.000 0 

47 . 000 

40.000 
4 Q .000 
50.000 





n - rw ' T ' w '^ TA»T 1 005 ^ , F1 1 M 005 1 , f 21 ( 1 ^ 05 ' , X ^ 1^01 ^ ^01 ,5 

Y5ir/t-pMr.rr)V 1 ( "^1 1 1 <' 1005' ,V2n''”1 M 00^ ^ ^ 1 025 

nTM^Mr; yr.fj ' 1 205 V, XI 1 1 M 005 1 1 1 1 M 005 ' ,X2?2 M 005 ^ 1 0 

0IM070 107 71 M 005 ^ , yp M no5 ^ , '"H'OTA 1 ^ 1 205 ' , 1 005 ) 

OAOA 7 DT '‘'’/1021 / 

1 X ’ l '^’^ OPA'^TOX '^XD ' 


I^TpTjm m 

Qpmptjm I X>3PP7 '-/OTOXT 1 i’ij'”)' 

Tf^Tpnm VP.I05'"1 

OXOPUT ' TXPO '^ 7 '^ IO-r ' 

I7PU? '•/■5IG5T2 

QTj^ t>up > r a-rlO L5705"' 1 ^ P? '' ' 

IX POP OAPLO^^'i 

QTJ7 pijP ' C A? LO L070 ’IT ? ( t''” 1 ' 

I7PTI? CAPLPX? 

OM 1 = ( WPTOH'^ 1 * . 45 " 6 ) / ^ 5? . 1 '74 * . 304 O ) 

S^';?=^4PT^7H7?*.4556^ /C52. 174 *. 50-1 5) 

5L1 =OAPLOX1 *.5015 
3L2 =0 AP L7. T-1 . 50 1 5 

0M1 1 =5L1 *5L1 *(5^M +072'^ 

5M1 2=c.p.-ip*7j,l *0T,o 

qwpl =0M -1 o 
am?2-'P'12*'^L?**2 

SKI 1 =5 2 . 1 * . 504 * ( OM 1 +0M 2 ) 

3K2?=5? . 174*.’^^5P.*OM?*OL2 

&AA= R-x-f f’PI 1 *"'M?2+P’.^22*0’'’1 1 )/(n^11 1 *‘^'’^22-22^ 2*2'“?1 ) 
xippl* 4*7 ( 70’/1 1 *2''P7+q’App*c;M'l \ i**2 1 /^‘^'♦’1 1 *2M22_5W1 2 

)**o_4*5K1 1 *0P?7/(nMl1*5M22-0’'''12*OM2i '))*’". 5 
W12Q = AAA-P'’5 
V/250=AAA+5P7 


'-/1 =f W1 7*^ )** . 5 
W?=(y 200 )**. 5 

1 61 , 20)'-.'1 , W 2 

>0 POPM A2 f 4 .. 6 1 V 

^070*1 =''0M1 1 +fP2'"1 1 -W 100*07.1 1 1*^721 +071 2 1 / (v/l 00*501 2 ) 
Xrqpi ■; lpi'^o*<^:M1 1 1**2 *nM2'’ /{'■!] 30*0'''1 2 '*>^2 )** .5 
' nox 0 ?=( 071 1 +( 571 1 - W 2 n 0 * 0 O 1 i 1 *( SM 21 + 5012I / fW 250 * 5 Oi 2 l 

. + 7o,9ri l_’-/230*roi 1 )**2*q^122/(W?50*0O1?1**2)**.5 

'oALCULA^P PIGPXVOOTORS 
PIG VI 1 =1 . 

PIGV1 2=1 . 

PIGV2 1' = (OKI 1 -Ml 00*PM 1 1 ) / C M1 50*''’'0 1 2 ) 

PI0V22=(5^1 1 -W'’0P*5M1 1 ) / ('42S0*r>wi 7 ) 

MPITP(l62,50lPT0V1 r,PIGV1 2,PTGVP1 ,PIGV2?. 

■7,0 POP''^ A'” f 4P1 4 . 6 ) 

CALOTILA7IOM O'^ PO^AL A1A7PIX 00’fAP0>IPX'!’0 MPPPPTCAL VALHPO 


TJ1 1 =1 .0/00X51 
TJ1 2 = 1 .0/00X5? 

TJ21 =fopi 1 -Miro^o^i 1 i/fMi 7 o*oMi 2*co:ri ) 

U 22 = ( 0^71 1 - M 77 n *+'"! 11 )/^ v ' 5 nn*qM 1 P '^’ OOX '^ P ' 


.non 
Ro , ono 

^5.000 

5?‘.ono 
S7 , ono 
5B.00O 

5Q.noo 
fo . 000 
G1 .000 
6? . 000 
6’^.QOO 
fi4.000 0 

r>5 . 000 
6fi.ooo 

67 . 000 

6'0 . ooo 

60.000 
70.000 

71 .000 

72.000 
7-5.000 

74.000 

75.000 

76.000 

77 . 000 

76.000 

70 . 000 

SO . 0^0 

51 .000 

62 . 000 
S? . 000 

54 . 000 

66 . 000 
86 . 000 

87 . 000 

36.000 r 

89.000 
00 . 000 

91 .000 

92 . 000 

93.000 

Q4.000 
Q5 . 000 
Q6.000 
Q7.0O0 

96.000 

90.000 
1 00 . 000 

101 .000 

102. 000 


ORSGI'lAl 53 61 

OF POOR QUAUTY 

40 A” I 4^^ 1 . 0 ' 

TMT>'i'^ DA '-’ 4 
00 6 c .1 = 1 , 1^01 

50 , 16'^'! "'TMn'f T ^ ^ vnno''’! r.T ^ , vDpoT'’{ n 

00 70 1 = 1 , 1001 

■^11 ( I ) = -(0'" 1 '‘'’Ll +0'"2 ’^"L1 '1 *-XDD0T1 ( T 1 
fr ( T 1 

70 CO’';'-' 7 7112 

T>0 6 -^ .T = 1 ,1001 

•P1 1 1 f,T 1 = D1 1 (,T )*'T1 1 +-‘^21 ' T 1*IT21 
r?pn2 (.J )-V'\ 1 rn«-ni '’+’^71 f.T 1*TI72 
0AL'^ULA'"T1 2"’A VALU-OO 17-007 <■ A O 
50 C07'’’TN"I7 

DO 60 .1 = 1 ,1701 

■'^U'll 1 (-T IsOINTWI *( ' 1 ^*7'1 1 1 1 ^V'l 

Y(.T)=-D1J71 1 r,Tl 
i=7I72(,T 1 
60 n 07 '^T 702 

CALL PC vy f X, V, 5,771 
00 140 J= 1 ,1001 

X 1 1 1 rjwxf.j'i 
511 1 

UO WPITLf 106, 1/.1 '''’'1 1 1 (.■T's ,51 1 1 f.T^ 

141 FOR'’ AO ( 22 14 . 61 
00 1 60 .T = 1 ,1001 

FLN1 1 ^ J 1=0 T7 ( T-7TM3(,T 1 1 ) *-^227 (j 1 /vio 

Yf.t1 = DIT71 1 (.Tl 
X(.T) = 7IMF(.T) 

160 CON'’’TM’JF 

GALL T'f’5Xy(X,Y,5,7PI'‘'1 
DO 146 J=1 ,1001 
X2??(J)=X(Jl 
5222 (J 1 = 5 (.11 

145 \VRITF(104, 1 a a; 1x 222 (J 1 ,5222fJ) 

146 ^ 0 RWAT( 2 FH . 6 ) 

CALC1JLA0F OYOTFO RF370^IDT1 f A1 1 

DO 210 1=1 ,1001 
XI (I1=L1 1 •’^Xl 1 1 (Il+ni2*572'5(l1 
X7fl)=Il21*51 1 1 rTl+U2?*5272(T) 

WRT'-' 2 ( 10 Q,? 001 -’’TMFU) f '^-1 (I^ ,X 2 (’Tl 
200 FORM AT (32 14 . 6) 

210 C02TTMU2 

DO 610 1 = 1 ,1001 

'"^FTAI (T) = 130.0m fl 1 / 6. 141 50265 
THR"A7f T ) = 1 60.n-*y2 (T 1 /" . 1 41 50265 

] 1 7^520l'‘'TMRfl 1 ,'"R2"'A1 (l1 ,'"R2''‘A7fl 1 
620 F0R9AT(6214.61 
610 Cn2'”TfFIP 

DO 71 1 = 1 , 1001 

F1 1 (1 1=- fov 1 *0L1 +0M2*0L1 I^'-XOOOTP ^ T 1 
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1 01 . Onn 
1 04 . 0^0 
1 05 . 000 
106.000 

1 07.000 
lOn.ono 

10Q.OOO 
110. OQO 

111 . 000 

112. 000 
Ul. 000 
114.000 

1 15.000 

1 16.000 

1 1 7.000 

1 18.000 
110. 000 
1 20.000 
1 21 .000 
1 22.000 
1 23.000 

124.000 

1 25.000 
126.000 

1 27.000 
1 2«.000 

1 29.000 

1 30 . 000 

1 31 .000 

1 32.000 

133.000 

134.000 

135.000 

1 36.000 

1 37.000 

1 38.000 

139.000 
1 40 . 000 
1 41 .000 

1 42.000 

1 43.000 
1 44 . 000 

145.000 

146.000 


i- 






OF POQil QUALll 


n-g 

15 J 


"’'71 ( J ^ = x '. vvnnriT ’ r * ( X 1 

71 oo'';'"i^ri8 

00 81 J=1 , 1orsi 

?1 1 1 (.1 )=T^i 1 r.T )-^ni 1 4-‘!i'7i (.1 

] f ,] Ml * riO 0 

81 On'MI'I”; 

8TA VAT/T8'’ Mol 

00 81 J=1,1001 

7'TIVI 1 fj U^X^lf'-ri *M-'PTwr./’.T 1 1 1 r.T ^ /vi 

Yf.jU-PTT^l 

X^Jl=?XMT^f.T) 

81 000TI\8I8 

OALL '^P37Y('X, 

00 143 J=1 , 10^1 
X1 1 1 (j)=x^n 
3111 (J) = 3^i) 

1 4"^ W2IT0^ 1 87 , 1 A2 IvM 1 t j 1 ^ 1 1 f ,1 1 

142 F0R’'*A'^('2'31/1 .61 
DO 1 «1 .T = 1 , 1001 

?U\ 1 1 ( ,T ) = 8I’'Tf '■fo*/' j 1 \ 1*-pppp ( j 'j /yo 

■^(J) = ‘^U011 fj} 

Xl'j) = '"X'^Ff'j) ’ 

181 COOT 1*108 

C A LL P3 77 ( X , V , 3 , HO I » 1 
00 14” -1 = 1 ,1801 
X222(-Tl=Xf.T'! 

3222(.I) = 3(J1 

1718 WPITXI' 111,1 48 1 Yopor.T > 3770 Cj ) 

14” ■^087A'‘'(7Xi4.‘f;i 

” ALOOLA'^Y r . Y '^ TT’l 988 P 0 H 88 XXt>ijti 
DO 21 1 1 = 1 ,1001 
XI ( I )=H1 1 -»-3l 1 1 (T 1+TI1 7*3070f 1 1 
X2 ( X ) = 021 -^3111 (I )+tiooif7/,no M 1 
WPXTEf 1 1 8, 308 MXfl8fl 1 XI M 1 "yo ( j ) 

300 T^08MAT(38U.61 ’ - 

211 C0HTTH1J8 

DO 400 1=1 , 1801 

"'HE''’A1 Cl 1 = 1 80. 8 !s-yi f T 1 /-^ . 1 /i 1 5qof,c, 

THFTACf T ) = 1 80, 8+X7 Cl 1 /- . 141 300A8 
WRITOf 1 1 3, 4,7n1'^T'':8r 1 1 ,'")F'"A1 Cl 1 a^Ct 1 
420 F0HXATC3X14.61 
400 OOHTI*-FIH 

pmr^p 

800 * 


appendix II 


COMPUTER PROGRAM FOR CALCULATING 
VALUES OF Hi FROM ANALYTICALLY 
INTEGRATED EQUATION 


ORIGii'IAf- 

OF POOR QUALlW 


1 .000 
2.000 
^.000 

4.000 

5.000 

6 . 000 

7.000 

8.000 
q.ono 

10.000 
11.000 
1 2.000 
15.000 
1 4 . 000 

15.000 

1 6 . 000 


^1=900 
T?=4 .■^4^'7 
)=0 
'•;i = . 

DO 50 1=1,400 

X 1= .0OOp*p*T (’T *(<^1 _x f 1 1 M /'‘fi *¥-0 

.+.OO0?’^^*'^I7fVM 1 1 1 /'.‘ri 

. -. 0002*i^*OI7O..'1 *ti ) /vn **5 
. - . 04 008 ^ W 1 * ^ 1 ^ T ) /V/1 **? 

. + . 04 *^-* OOOfW1 *'^1 i /V1 **'’ 

50 TCT+1 )=?('!'> + . 

DO 10 T=1 ,400 

1 0 V/DI7'J=' h 0,0 , 1 00 ■' ^ I ^ , PTf ( 1 1 

1 00 AO f OD 1 4 . 

O'^O? 
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